Quantum Evaporation of a Bose-Einstein Condensate 
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We show that a Bose-Einstein condensate emits atoms, if either the condensate wave function, 
or the scattering length of the atoms depends strongly on time. Moreover, the emission process 
is coherent and atoms can oscillate back and forth between the condensate and the excited states. 
Inspired by recent experimental results, we present results of simulations of the response of a Bose- 
Einstein condensate to a very rapid change in the scattering length. The possibility of molecule 
formation is also discussed. 
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Introduction. — There are several examples where the 
coupling of a classical field to a quantum field leads to the 
production of quanta of the latter. One of the best known 
examples is Hawking radiation, where the curved space- 
time provides the classical 'background' that couples to 
a quantum field. Another example is bremsstrahlung, 
where an electron scatters with a nucleus, and emits 
a quantum of the electromagnetic field, i.e., a photon. 
From a quantum mechanical point of view, the emission 
of radiation by an antenna is also described by the cou- 
pling between a classical current, and the photon field. 

An example of a classical, albeit complex field, which 
is subject to a high level of experimental control, is the 
macroscopic wave function of a Bose-Einstein conden- 
sate. Recently, this experimental control could even be 
extended to the interactions between the atoms |[]-[|, 
which are at the temperatures and densities of interest 
solely determined by the s-wave scattering length. Mak- 
ing use of a so-called Feshbach resonance , one is able 
to vary the scattering length of the atoms to any pos- 
sible value, by tuning the magnetic bias field. This ex- 
perimental degree of freedom was recently exploited to 
study the condensate collapse Jf|, which was first ob- 
served in the experiment by Bradley et at Q. Roughly 
speaking, such a collapse occurs if the condensate con- 
tains so many atoms that the mean-field interactions ex- 
ceed the energy-level splitting of the external trapping 
potential . Clearly, the wave function of a collapsing 
condensate is an example of a time-dependent classical 
field, since it undergoes very rapid and violent dynamics 
ill- 

Another type of time dependence arises if the scatter- 
ing length changes on a time scale that is fast compared 
to the collective modes of the condensate, i.e., fast com- 
pared to the inverse of the frequencies of the external 
trapping potential. 

It is the main purpose of this Letter to point out that 
both kind of time dependences, can cause the transfer of 
condensate atoms to excited states. Theoretically, this is 
described by the coupling between a classical field, i.e., 
the condensate wave function, and a quantum field which 



describes the atoms in the excited states, similar to the 
examples mentioned above. Moreover, due to the fact 
that this quantum evaporation process also contains a 
coherent part, the atoms can oscillate back and forth be- 
tween the condensate and the excited states, in a kind of 
multimode Rabi oscillation. This condensate loss mecha- 
nism has therefore the peculiar feature that at later times 
less atoms are lost from the condensate. This behavior is 
completely different from conventional loss mechanisms, 
such as dipolar relaxation and three-body recombination, 
which are characterized by a rate coefficient and there- 
fore always lead to more atom loss after a longer time 
evolution. In a very recent experiment by Claussen et 
al. |l4|] , the number of condensate atoms was measured 
after a very rapid change in the scattering length. In this 
experiment, it was found that the number of condensate 
atoms indeed increases with time in some regimes, and 
therefore we believe that our theory might offer an ex- 
planation for these experiments. 

Quantum evaporation. — A convenient starting point 
for our discussion is the second-quantized hamiltonian 
for the system. It is given in terms of the Heisenberg 
creation and annihilation operators, denoted by if) ''(x, i) 
and tp(pc,t), that create and annihilate an atom at posi- 
tion x and time t, respectively, and obey the usual Bosc 
commutation relations at equal time. The hamiltonian 
reads 

dx ^(x, £)«oV>(x, t) 



H 



+ ^ Ll ydx V t (x,i)^ t (x,t)^(x,i)V(x,t). (1) 

The single-particle hamiltonian Tlo contains the kinetic 
energy of the atoms and the external trapping potential. 
Its eigenstates and eigenvalues are denoted by %n (x) and 
e n , respectively. The strength and sign of the interac- 
tions are fully determined by the two-body T(ransition) 
matrix, T 2B (t) — ATra(t)h 2 /m, where m is the mass of 
one atom. Note that the s-wave scattering length a(t) 
is explicitly allowed to depend on time, which is exper- 
imentally realized by tuning the magnetic field near the 
Feshbach resonance. 
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We split the operators in a condensate part, denoted 
by (j>(x, t), which will be treated as a complex classical 
field, and a part that describes the fluctuations, denoted 
by V>'(x, t). We insert the separation ip = 4> + ip' into the 
hamiltonian in Eq. (|l|), and keep terms up to quadratic 
order in the fluctuations. This is known as the Bogoli- 
ubov approximation, and means physically that we as- 
sume that most of the atoms are in the condensate. This 
assumption is reasonable for the zero-temperature appli- 
cations under consideration here. The hamiltonian now 
consists of three parts, Hqp + Hb + flint j where Hqp 
is the usual Gross-Pitaevskii energy functional, and Hb 
contains the terms that are quadratic in the fluctuations. 
The interaction between the condensate and the quantum 
fluctuations is described by fl; n t , which, in first instance, 
is given by the sum of 



dxcj>*(x,t) Ho +T 2B (i)|0(x,i)| 2 ip'(x,t) 



and the hermitian conjugate expression. 

Physically, we want the classical field (f>(x,t) to de- 
scribe the low-energy part of the system, i.e., the low- 
energy single-particle states, whereas the fluctuations are 
the high-energy excited states. Depending on the physics 
of the specific application, we introduce a cut-off between 
these two parts of the system. This implies that we drop 
the terms containing Ho in if; n t, which now takes the 
form of a coupling between the classical 'current density' 
J(x,t) = T 2B (i)|0(x,i)| 2 0(x,i), and the quantum field 
ip'(x, t). Because of this coupling, a perturbation of the 
classical field can result in the production of quanta of 
the quantum field, i.e., atoms can be transferred from the 
condensate to the excited states. To study this process, 
and to derive a rate equation for the number of atoms in 
the condensate, we have to solve the Heisenberg equation 
of motion for the quantum field operators, given by 



ih 



d_ 

dt 



H -2T 2B (t)\cj ) (x,t)\ 2 ^'(x,t) = J(x,t), (2) 



where the hermitian conjugate expression holds for the 
creation operator ^''(x, t). In this equation of motion, 
we have neglected the anomalous parts of the hamilto- 
nian Hb, since these describe the collective motion of 
the condensate, and are supposed to be included in the 
condensate wave function cf>(x, t). 

The Heisenberg equation of motion in Eq. (J2J) is most 
easily solved by introducing the retarded Green's func- 
tion G( + )(x, <; x', t') by means of 

ih^--H - 2r 2B (t)|0(x,i)| 2 lG(+)(x,t;x',t') 
ot J 

= h8(x-x')5(t-t'), (3) 

with the boundary condition (x, i; x', t') — for 
t < t' . Physically, this Green's function describes the 



propagation of the atoms in the excited states, in the 
absence of the interaction H mi . Therefore, we have that 

iG(+)(x, t; x', t') = 9(t - t')([i>'(x, t), ^'t(x', t')]) J=0 . (4) 

Because of the coupling with the classical 'current den- 
sity' J(x, t), the operator for the fluctuations acquires a 
nonzero expectation value, given by 

(ip'(x,t)) = ~ Jdt' Jdx' G (+) (x,i;x',i')J(x',0- (5) 

Assuming that initially all the atoms are in the conden- 
sate the density of the noncondensed atoms is given by 
n'(x, t) = I (ip'(x, t)) | 2 , and the rate equation for the num- 
ber of atoms N c in the condensate now reads 



dN c (t) 
dt 



d f 

dtj y ' 



— rr dx dt' dx'Jm 



T 2 «(tMx,i)|V(x,i) 



GW ( Xj t ; x', t')T 2B (t>)\4>(x>, i')| 2 0(x', t')] . (6) 

This equation is our most important result, and de- 
scribes the change in the number of condensate atoms due 
to a strong time dependence of the condensate wave func- 
tion, and/or the scattering length of the atoms. Since we 
are treating the quantum field ip 1 as noninteracting, the 
process is coherent. In particular, oscillations of atoms 
between the condensate and the thermal cloud can oc- 
cur. Clearly, the rate equation is nonmarkovian, and on 
short time scales there is no conservation of energy, due 
to the Heisenberg uncertainty principle. However, in the 
markovian limit, where the energy of the ejected atoms is 
taken much larger than the energy of a condensate atom, 
the rate equation takes the form of Fermi's Golden Rule 
for the process of an elastic collision between two conden- 
sate atoms, where one atom is ejected from the conden- 
sate and one atom is stimulated back into the condensate 
|]l3| . In particular, this means that in equilibrium there 
is no correction to the usual Gross-Pitaevskii equation. 
Presently, we are mainly interested in the new features 
that arise due to the nonmarkovian nature of the process 
under consideration. 

Multimode Rabi oscillations — As an example, we dis- 
cuss the response of the system to a change in the scat- 
tering length that is too fast for the condensate to react 
dynamically. More specifically, we do our calculations for 
the experimental parameters of Claussen et al. fl4|| . In 
these experiments, 85 Rb atoms are confined in a cigar- 
shaped trap with radial frequency uj r /2i: — 17 A Hz and 
axial frequency lo z /2tt = 6.8 Hz, and the Feshbach res- 
onance at 156.9 (G)auss is used to vary the scattering 
length very rapidly. This is achieved by a trapezoidal 
pulse in the magnetic field, which means that the field 
is ramped linearly to a certain value in a time t r i se , then 
held for a time thoid, before ramping back to the ini- 
tial value. The initial and final values of the magnetic 
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field are chosen such that the scattering length is ini- 
tially equal to zero, and is large and positive during the 
hold. The rise time and the hold time are typically of 
the order of microseconds, and therefore the shape of 
the condensate wave function hardly changes during the 
pulse, but remains the ground state of the trap, i.e., a 
gaussian. However, the phase of the condensate wave 
function changes considerably, and is in a gaussian ap- 
proximation given by |L3| 

ft M j. m dq i^ 
9 ^) + ^L. qj{t ) dt ' 

with qj(t) the width of the gaussian in the three spa- 
tial directions. Including both global and local phases of 
the condensate wave function in our calculations is im- 
portant, since, roughly speaking, these phases determine 
the energy of a condensate atom. The variational param- 
eters qj (t) turn out to obey Newton's equations of motion 
H . The equation of motion for the global phase 6*o(i) is 
determined by the condensate energy, which includes the 
effects of the evaporation process. 

For the retarded propagator of the ejected atoms, we 
use the expansion 

G(+)(x,t;x',t') = -^-0 V-n(M') 

xx^x^xOe-^— *'>, (7) 

in which the sum is over all eigenstates, except for 
the ground state, which is already contained in 0(x, t). 
The coupled equations for the expansion coefficients 
a n , m (t, t') are found from the equation of motion for the 
Green's function in Eq. (||). These coefficients clearly 
obey the initial condition a njm (i,t) = <5 n!m at the initial 
time, since all the atoms are then in the condensate. 
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FIG. 1. The fraction of atoms in the condensate as a func- 
tion of the hold time, for a fixed rise time of t T i ac = 12.5/is. 
The solid line corresponds to initially N c (0) — 6100 atoms 
in the condensate, whereas the dashed line corresponds to 
N c (0) = 16500 atoms. The scattering length during the hold 
is equal to a = 2000ao . 
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FIG. 2. The fraction of atoms in the condensate, as a func- 
tion of the rise time, for different hold times, for a condensate 
of initially N c (0) — 16500 atoms. The solid line corresponds 
to thoid = 1/is, the dashed line to ihoid = 5/xs, and the dotted 
line to thoid = 15/is. The scattering length during the hold is 
equal to a — 2000ao . 



In Fig. [l] we show the fraction of atoms in the con- 
densate as a function of the hold time, for a fixed rise 
time of i r iso = 12.5/is. The initial number of conden- 
sate atoms is N c (0) = 16500 atoms, and A c (0) = 6100 
atoms, respectively. The pulse is such that the scattering 
length is equal to a = 2000ao during hold. Here, a is 
the Bohr radius. A significant fraction of the atoms is 
transferred to the excited states in both cases, and part 
of this fraction can come back into the condensate after 
some time, as seen from Fig. [IJ. The atoms come back 
faster in the case of the largest initial number of atoms, 
which is caused by the fact that the coupling between the 
condensate and the excited states is proportional to the 
number of atoms. Note also that there are various fre- 
quencies in the curve, which displays the fact that we are 
dealing with oscillations of atoms between the conden- 
sate and several excited states. In Fig. the fraction of 
atoms in the condensate as a function of the rise time is 
displayed, for various hold times. Clearly, the number of 
atoms increases with the rise time, over some range. This 
can, of course, not be understood from the viewpoint of 
a loss mechanism characterized by a rate constant, such 
as three-body recombination, or dipolar decay. 

Discussion. — Comparing Fig. ^ to the experimental 
results of Claussen et al. [[l4| , we notice that experimen- 
tally the number of atoms always decreases with increas- 
ing hold time. An explanation for this behavior is the 
presence of three-body recombination. The comparison 
of Fig. |^ with the experimental results leads to the same 
conclusion. The minima in this figure are seen to shift to 
the left, as in the experiments, and also occur at the cor- 
rect value of £ r isc- However, experimentally they become 
lower with increasing hold time, 
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FIG. 3. The fraction of atoms in the condensate as a func- 
tion of the hold time, for a fixed rise time of t r i ae = 12.5/j.s. 
The solid line corresponds to a — 2000ao and the dashed line 
to a = 4.9 x 10 5 ao during hold. The initial number of con- 
densate atoms is equal to N c (0) = 16500 atoms. 

which can also be explained by a background loss mech- 
anism independent of the loss mechanism discussed here. 
Such an additional loss mechanism leads also to the 
smoothing of the higher-order oscillations seen in our 
numerical results, which are not observed in the exper- 
iments. Unfortunately, near a Feshbach resonance the 
behavior of the three-body recombination is not known 
as a function of the magnetic field. Therefore, it turns 
out that a more detailed comparison with experiment is 
impossible at this point. 

A popular model for the description of a Feshbach reso- 
nance involves the coupling of an atomic field to a molec- 
ular field M,E6|. In such a model Rabi oscillations can 
also occur. To the extend that we may neglect quantum 
evaporation, we can solve this model essentially exactly 
using the framework developed here, by making the sub- 
stitution T 2B \(f)\ 2 (l)*i>' -> g(4>*) 2 '4> m in H mt , where ip m 
annihilates a molecule. Note that we do not have to take 
into account the anomalous averages of the atomic oper- 
ators, since these are effectively already incorporated by 
using the experimental values of the coupling constants. 
Doing so would lead to a double counting of the inter- 
action effects Jn|. In the Heisenberg equation of motion 
for the molecular field operator, we have to incorporate 
the detuning from the Feshbach resonance, defined as 
the energy difference between the energy of the bound 
molecular state with respect to the atoms. It is given 
by e(t) = n{B(t) - B )/AB, where B and AB denote 
the position and the width of the resonance, respectively. 
The prefactor k is fixed with the knowledge of the Zee- 
man effect of the molecule with respect to the atoms 
p8[ , and the coupling g is subsequently determined to 
reproduce the correct strength of the Feshbach resonance 
p5[ . Since the 85 Rb molecular bound state responsible 
for the Feshbach resonance can not be trapped magnet- 
ically, we perform our calculations for an optical trap, 



with the same frequencies as used in our previous calcu- 
lations. In particular, this means that the energy levels 
of the atoms and molecules are the same. We perform 
our calculations for different hold times, but a fixed rise 
time of t r ise = 12.5/is. In Fig. || we present the result 
of our calculations, for two different values of the mag- 
netic field during hold and for a condensate of initially 
N c (0) = 16500 atoms. Note that for the experimentally 
relevant case of a = 2000ao essentially no atoms are con- 
verted to molecules. This can be understood from the 
fact that for this value the detuning e is much larger 
than the effective coupling, i.e., e 3> g\fn, with n the 
density of the atomic condensate. This means that the 
amplitude of the multimode Rabi oscillations between 
the atomic condensate and the molecular gas, which is 
of order O ((g-^/n/e) 2 ) , is very small. In order to con- 
vert a significant fraction of atoms to molecules, we have 
to consider magnetic fields very close to the resonance. 
In Fig. H we also show the results of our calculations 
for a = 4.9 x 10 5 ao, which corresponds to a magnetic 
field that is 0.01 G above resonance. Due to the collec- 
tive motion of the condensate the coupling between the 
atomic and the molecular trap states decreases with time, 
which leads to a damping and dephasing of the oscilla- 
tions. Since we are dealing with an enormous effective 
scattering length, the decay through three-body recom- 
bination, and spin-flip processes becomes important. In 
our calculations, we have however, for the same reasons 
as before, not taken these decay processes into account. 
They will lead to a further dephasing of the Rabi oscil- 
lations. 

In a future publication we intend to study the two 
pulse experiments conducted recently and look at 
the properties of the ejected atoms in more detail. 

It is a pleasure to thank Servaas Kokkelmans, Neil 
Claussen, and Eric Cornell for helpful remarks. 
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